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Abstract. We have searched for periodicity in the photo- 
metric and spectroscopic variations of the enigmatic star 
HD 101584 and found a long-term variability on a typical 
time scale of 1700 days (4.7 years) and a highly significant 
period of 218 ± 0.7 day. This period is most prominently 
present in the photometric indices which are a measure 
for the Balmer jump (Geneva d and Stromgren c?). The 
Doppler velocities of the high-excitation photospheric ab- 
sorption lines (He I and C II) seem to be variable with the 
photometric period. Our data favor the 218 day period 
for the Doppler velocities with a small probability that 
the true period is 436 days. We argue that HD 101584 
is a close (highly) eccentric 218 day binary system with a 
low-mass unseen secondary. The photometric and Doppler 
variations are attributed to changes in the velocity law and 
mass-loss rate of the stellar wind which lead to asymmet- 
ric line profiles and a phase dependent Balmer disconti- 
nuity. Binary interaction is responsible for the changes in 
velocity law and mass-loss rate leading to the observed 
phenomena. 

Key words: techniques: photometry - techniques: radial 
velocities - binaries: close - stars: HD 101584 - stars: mass- 
loss - stars: AGB and post-AGB 



1. Introduction 



Since Humpreys ( 1976| ), on the basis of limited data, sug- 
gested that the peculiar supergiant HD 101584 is a binary 
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system with an orbital period of 3.5 years, it has been an 
open question whether the binary nature could be proven 
independently. The last two decades no convincing im- 
provements have been made in detecting a periodicity. 
Here we will show that the optical photometry shows a 
highly significant periodicity of 218.0 ± 0.7 days and that 
the same period is probably present in the Doppler veloc- 
ities of high-excitation absorption lines. 



Parthasarathy and Pottasch ( 1986| ) and Trams et 
al. (1991) used the detection of infrared excess in the 



IRAS bands as an argument for the evolved nature of 
HD 101584. They interpreted this as radiation from the 
AGB ejecta and argue that HD 101584 is in the post- 
AGB phase. The galactic latitude of 6° is too low to decide 
between population I and II, but the observed radial ve- 
locity of HD 101584, derived from the CO millim eter line 
emission 

(w,co = 50 km s-i , Trams et al. \l99(\ ) cannot 
be fitted with the Galactic rotation curve if the star is a 
massive population I object. This strongly points to the 
population II nature of the object. 



Bakker et al. (1996) have made an extensive study on 
the nature of HD 101584 by combining multi-wavelength 
photometry and spectroscopy and found that the system 
contains a late B-type post-AGB stars which experiences 
a very high mass-loss rate (M ~ 10^* M0 yr^^). Due to 
the high density, the wind has a low degree of ionization, 
close to the star, the ionization conditions in the wind are 
similar to those in the photosphere of a F-type supergiant. 
This is observed as a low ionization shell spectrum. Due 
to the numerous UV spectral lines in the shell spectrum 
the UV continuum is not observed but obscured by a cur- 
tain of absorption lines (Bakker et al. 1994). In the red 



spectrum the shell spectrum produces less spectral lines 
and the late B-type star is observed by the presence of 
high-excitation absorption lines (e.g., He I and CII). 
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The most obvious way to look for binarity is to study 
the Doppler velocities of spectral lines, which gives the 
orbital period and the velocity amplitude of the ob- 
served star. This method cannot easily be applied to 
HD 101584 because of the velocity-stratification in the 
pseudo-photosphere. The Doppler velocities of absorption 
lines are not only dependent on the orbital phase, but also 
on the wind structure, line strength and central wave- 
length of the line (Bakker |l993| ; Bakker et al. |1996D . A 



We have calibrated the Stromgren photometry of 
HD 101584 using the method described by Spoon et al. 



radial velocity curve on the basis of high-excitation pho- 
tospheric absorption lines may prove the binarity. 

We have looked for periodicity in the photometric data 
of HD 101584. One of us, CW, has observed HD 101584 
over many years, and obtained 89 observations in the 
Geneva system. The second, independent, dataset we have 
used is from the Long-Term Photometry of Variables 
(LTPV) project. These 75 observations have been ob- 
tained in the Stromgren system, were the y-band has been 
transformed to the Johnson l/-band (Vj). A description 
of the observations is given in Sect. ^. We have combined 
the Geneva F-band with the LTPV y-band to obtain 157 
Johnson Vj magnitudes over 3337 days (4.1 years). Using 
the Phase Dispersion Minimization method (Sect. ||), we 
have detected a 218.0 ±0.7 day periodicity in Vj (Sect. ^. 
In Sect. ^ we will show that the 218 day period is likely 
present in the Doppler velocities of high-excitation photo- 
spheric absorption lines, and discuss the results in Sect. ^ 

2. The observations 

2.1. Photometry in Geneva seven color system 

CW has collected photometry in the Geneva seven color 
system ([/, Bl, B, B2, VI, V, G) spanning 1511 days 
(4.1 years). These observation were obtained with the 
70cm Swiss telescope at ESO, La Silla, Chile. The error 
estimate on the V magnitude is 0.003 mag. A complete de- 
scription of the Geneva seven color system can be found in 
Rufener and Nicolet (1988). We will use the Geneva pho- 
tometric indices B2 — V\, d and g which are a measure of 
the effective temperature, the Balmer discontinuity, and 
the blocking by absorption lines respectively. 

2.2. Photometry in Stromgren uvhy system 

We have obtained data from the Long-Term Photome- 



try of Variables (LTPV) project (Sterken 1983) covering 
10.6 years. All observations were taken using the Dan- 
ish 50cm telescope on La Silla, and are made in the 
Stromgren uvby system. The photometry of the y-band 
has been transformed to the Johnson Vj-band usi ng th e 
relation between those two bands (Manfroid et al. 1991 ). 
The data has b een pu blished by Manfroid et al. (1991) and 
Sterken et al. ( 1993 ). We will use the Stromgren indices 
(6 — y)o, c\, and mj, which are a measure of the effective 
temperature, Balmer jump, and blocking by absorption 
lines respectively. 



(1994). The LTPV program makes use of comparison stars 
to increase the accuracy of the data. Just before and just 
after observing the program star, comparison stars are 
observed. Comparison star A is HD 102113 (AOIV) and 
comparison star B is HD 101735 (G5III). By assuming 
that the two comparison stars are intrinsically constant in 
magnitude we can calculate a correction term by which the 
observed magnitude of the comparison stars differs from 
their average value over the whole dataset. The observed 
magnitude of the program stars is corrected with this cor- 
rection term to obtain an "atmospheric extinction" free 
intrinsic magnitude (Eq. |l]). In order to correct the data, 
the comparison stars have to be observed within 30 min- 
utes before or after the program star. Observations of the 
program star within three minutes of each other were av- 
eraged to one magnitude. 



_ m,A + mB , < mA + rriB > ^ 

^corrected — ^observed i 



where < rriA + ms > is the mean magnitude of stars A 
and B averaged over all observations. 

2.3. The Johnson Vj-band 

The filter response function of the Geneva F-band and 
the Stromgren y-band, as converted to the Johnson Vj- 
band, are almost identical. For this work we have made 
use of this and compiled a dataset of Johnson Vj-band 
photometry containing the Geneva y-band and the cali- 
brated Stromgren y-band. This dataset contains 157 data 
points and will be referred to as the Johnson Vj-band 
dataset throughout this article. 



2.4- Doppler velocities 

High-resolution CAT/CES spectra (Table 0, 17 spectra) 
have been obtained by CW in the wavelength region of 
He I (5876 A) and C II (6579 and 6583 A). Although there 
are many lines in the optical spectrum of HD 101584, the 



majority is formed in the stellar wind (Bakker et al. 1996) 
and only the high-excitation absorption lines (e.g., He I 
and C II) are from the photosphere of the underlying B- 
type star. Here we will try to correlate the photometric 
period with the Doppler velocities of these high-excitation 
absorption lines. All He I lines were fitted with a Gaussian 
to derive the central velocity. The C II lines are too weak 
for a reliable Gaussian fit and the central velocity was de- 
termined with an eyeball fit and the equivalent width by 
integrating over the line profile. If more than one velocity 
is available at a given date, the velocity used for our anal- 
ysis was based on the quality of the spectra (Table . A 
typical error in the velocity is 3 km s^^ . All velocities are 
hclioccntrically corrected. 
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Table 1. Heliocentric radial velocities of the C 11 (6579 and 6583 A) and Hel (5876 A) absorption lines. 



Date 


JD (244+) 


Id. 


^obs 


W 




dd/mm/yy 






[A] 


[mA] 


[km 


17/02/89 


7575.85 


Cll 


6578.69 


35 


48 






Cll 


6583.55 


25 


48 


23/01/90 


7915.43 


Cll 


6578.87 


38 


57 






Cll 


6583.45 


14 


46 


02/07/90 


8074.75 


Hel 


5876.98 


111 


52 


14/03/91 


8329.65 


Cll 


6579.01 


30 


55 






Cll 


6583.61 


24 


45 


29/01/92 


8650.86 


Cll 


6578.61 


? 


45 






Cll 


6583.52 


19 


49 


17/04/92 


8729.50 


Cll 


6579.18 


44 


53 






Cll 


6583.92 


26 


49 


18/01/93 


9006.0 


Hel 


5876.32 


158 


55 


18/01/93 


9006.0 


Cll 


6578.81 


48 


55 






Cll 


6583.52 


53 


50 


11/02/93 


9030.0 


Cll 


6579.05 


86 


64 






Cll 


6583.77 


45 


59 


12/02/93 


9031.0 


Hel 


5876.49 


143 


62 


22/04/93 


9089.69 


Hel 


5876.81 


155 


59 


22/04/93 


9098.95 


Cll 


6579.30 


48 


56 






Cll 


6584.12 


29 


56 


14/05/93 


9121.54 


Cll 


6579.52 


33 


60 






Cll 


6584.27 


31 


56 


15/05/93 


9123.46 


Hel 


5876.90 


162 


57 


19/07/93 


9188.46 


Hel 


5877.10 


145 


57 


20/07/93 


9189.50 


Cll 


6579.69 


43 


57 






Cll 


6584.55 


35 


59 


14/12/93 


9335.83 


Cll 


6578.89 


31 


55 






Cll 


6583.61 


22 


51 



Remark 



very noisy 

very noisy 
very noisy 
very noisy 
noisy 



nice profile 
nice profile 
nice profile 
polluted 



V.JD 

[km s^""" 



48.0 



57.3 



51.8 



52.6 



55.3 
52.4 



64.1 



61.8 



asymmetric 59.2 
56.3 



58.5 



57.3 



57.0 
57.9 



55.3 



3. Phase Dispersion Minimization method 

Searching for periodicity with Fourier analysis techniques 
one looses information in trying to fit a sinus shaped curve 
to the mean light curve. This loss of information does not 
occur using the Phase Dispersion Minimization method 
(PDM). The PDM method was first proposed by Lafler 



and Kinman (1965) for finding the period of a RR-Lyrae 
variable. Stellingwerf (1978) made a generalization of this 
PDM method and applied it to the double-mode Cepheid 
BK Centauri. This method looks for periodicity in the 
light curve by means of an auto-correlation method that 
does not require the light curve to have any pre-known 
shape, but only that it repeats itself every period. 

A set of N photometric observations contains N pairs 
of Julian dates JD and observed magnitudes: {JDi,Vi : 
i = I, N). Let the variance of the magnitude, cr^, be given 

by 



N-l 

where V is the average magnitude V ^ J2 Vi/N. 



(2) 



The dataset is folded over a given period 11, to 
determine the phase of the observation dates (f> = 
{JD — J Do) /n mod 1. Dividing the phase in Nb bins and 
covering every bin Nc time with a bin-offset of l/{NbNc) 
we have M — Ni,Nc dependent samples of the mean light 
curve. 

For each j-th sample containing M points the variance 
Sj (Eq. H) can be computed analog to Eq. |^. 



Uj - 1 



(3) 



The overall variance for the given period 11 is then given 

by 

- ^^^-^^ (4) 
" E-f in, - 1) ^ ^ 

we use the ratio 

e^i (5) 

If the phase diagram shows a correlation for a given period 
n, the overall variance will be smaller than the dataset 
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variance a^. This gives a minimum in the 0-statistics 
(Eq. H). If the phase diagram does not show any corre- 
lation but only a random distribution of magnitudes, the 
overall variance will equal the dataset variance , and 
8 will be close to 1. 

The strength of the PDM method is that no pre-known 
shape of the light curve has to be given. For the photomet- 
ric data we taken {Ni,,Nc) = (10,2). This means that we 
split the light curve in 10 bins, each covering 0.1 part of 
the photometric phase. Each data point will be used twice. 
The 8-statistics was calculated for 0.0 to 0.01 days in 
frequency steps of 1 x 10^^ days^^. 

4. Photometry 

The long-term photometric variability with a time 
scale of 4-7 years 

The photometric data of HD 101584 shows a long-term 
variation on a typical time scale of 1700 days (4.7 years). 
Fig. |l| shows the U, B, and V broad band filters of the 
Geneva dataset (until JD = 2448349.6). While the U- 
band has a minimum near JD — 7700, the y-band has a 
maximum. In the next sections we will discuss the 218 day 
photometric period and find a similar behavior which is at- 
tributed to changes in the wind structure and successively 
in the Balmer jump. It is suggestive to think that this vari- 
ability is periodic, but our dataset is to limited to reach 
a firm conclusion about that. We note that the Kelvin- 
Helmholtz time scale for a post-AGB star is of the order 
of several decades. This might indicate that the origin of 
the variability has to do with the release of gravitational 
energy of a part (envelope ?) of the star itself. However 
the primary goal of this study it to look for photomet- 
ric variations which might indicate that HD 101584 is an 
interacting binary system. It is unlikely that the 4.7 year 
variability is due to binary motion as in such a wide bi- 
nary system the mutual interaction will be very small. Our 
main goal is to understand the large scatter around this 
long-term variation: the 218 days period. 

4-.2. The 218.0 ± 0.7 day periodicity in the Johnson Vj- 
band 

In order to look for a short-time periodicity on a long-term 
variability we should first try to correct for the long-term 
variability. However the shape of the intrinsic long-term 
variation is unknown, and we have tried to correct for this 
long-term variation without satisfactory results. The rea- 
son for this is that the PDM method is very suitable for 
finding a short-term variation on top of a long-term vari- 
ation. By correcting the photometry we would pollute the 
light-curve with and "arbitrarily" long-term light curve 
without gaining much in sensitivity of detection a short 
time periodicity. 

Using the PDM method as described in the preced- 
ing section we find a H = 218.0 ± 0.7 day period in the 



Johnson Fj-band of HD 101584 (Fig. |). The e-spectrum 
does not only show the primary frequency, /o — 456.446 x 
10~^ days~^, but also at least five sub-harmonics and two 
beats with the year. This is a highly significant period in 
the Johnson Vj-band. Table ^ gives a list of the identi- 
fied features in the 8-spectrum. Using this period we can 
make a least-square-fit to the date using a sine function 
(Table I). 

The many unidentified features in the 8-spectrum are 
likely due to aliases as a result of the irregular spaced grid 
of observation dates. We have made a synthetic light curve 
using a sine-function and the photometric parameters of 
Table ^. Such a synthetic light curve on the same grid of 
observation dates should show the same features in the 
8-spectrum. If this is not the case then this could mean 
that there is another period present in the data. Fig. ^.b 
shows the 8-statistics of the best-fitted-sine-function (Ta- 
ble H) on the grid of observation dates. All features are 
present but the power in the 8-spectrum is much more 
concentrated in the (sub-) harmonics of the sine function. 
This tells use that it should be possible to obtain better 
results using another shape of the synthetic light curve 
(e.g., sawtooth), or subtracting the long-term variability. 

Finally we have computed a synthetic spectrum on 
a daily grid between the first and last day of observa- 
tions present in the data set. The resulting 8-spectrum 
(Fig. ^c) shows very strong features at the primary fre- 
quency and at least seven sub-harmonics, no harmonics. 
The beats with the year are no longer present. The peak at 
primary frequency shows lobes on both sides correspond- 
ing to a frequency of the total time span of observations. 

A 218 day period in the 8-spectrum does not unam- 
biguously means that the true period is 218 days. In the 
case the photometric periodicity is due to the rotating ge- 
ometry of a binary system the true period can be twice the 
photometric period. As long as it is not clear what causes 
this periodicity it would be premature to choose between 
a 218 and 436 day period. 

Fig. ^ shows the Johnson V magnitude as function of 
phase (based on 218 and 436 day period) and Julian 
date. By looking at the light curve as function of Julian 
date we see that the magnitude follows the sine function 
quite well. We note that there are strong deviations from 
the mean light curve, particularly at its maximum. Un- 
fortunately the photometric light curve does not signifi- 
cantly show that two successive 218 day periods behave 
differently. The question whether the true period is 218 or 
436 days remains therefore open. 

4-3. The 218.0 ±0.7 day period in the photometric indices 

In order to find the origin of the photometric period we 
searched for periodicity in the other photometric bands 
and photometric indices. As the dataset is about half the 
size of the Vj-band dataset the PDM method does not 
give any useful results and we folded the photometry over 
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Fig. 1. The magnitudes of the Geneva broad band filters as function of Julian date show that there is a long-term variation in 
the photometry on a typical time scale of 1700 days (4.7 years). 



Table 2. Identification of the features in the O-statistics of 
the Johnson Vj-band. 

u Identification Period 11 

[IQ-^^daysJ^] [days] 



78.099 




213.40 


91.012 


fo/5 


219.75 


114.783 


/o/4 


217.80 


152.080 


/o/3 


219.18 


181.525 


fo - fs 


219.63 


229.827 


fo/2 


217.55 


456.446 


fo 


219.08 


735.732 


fo + fs 


216.47 


average 


fo 


218.0 ±0.7 


fs is the the 


sidereal year of 365.256 days 


J. Parameters of best fitted 


sine function. 


Parameter 


Best 


itted value 


Period 11 


218.0 


± 0.7 days 




6.934 


±0.008 mag 


Amplitude 


Al 0.016 


± 0.008 mag 


JD{(j>phot - 


= 0) 2445044.576 



i'phot 



is defined where Vj = Vj and dVj/dt < 



a 218 day period and 436 days. In the case of 436 days we 
find that the two successive sines (of each 218 days) do not 
differ significantly. Which means that our data favor a true 
period of 218 days. In order to increase the accuracy of the 
light curve we will fold the data with a 218 day period. No 
clear phase dependence is observed in any of the separate 
photometric bands (with the exception of the Johnson Vj- 
band because it contains much more data points). This is 
presumably due to the small amplitude and the irregu- 
lar behavior of the light curve. To show the presence of 
photometric variations we will first look at the standard 
deviation (a) of the photometry for each band and index. 

Table ^ gives for each photometric band (in order of 
increasing wavelength) the average magnitude (M) and 
the standard deviation (ctm)- The last column gives the 
standard deviation of that band or index for the compar- 
ison stars (cta/b) and is used as a measure of the error in 
the dataset. Here we note that the standard deviation on 
the average magnitude of a photometric band or color is 
significantly larger than the error. Fig. ^ shows that the 
standard deviation increases to the blue. Near the y-band 
the variations are on a 3(T level, while due to an increase 
of the error the variations decrease to the blue. From this 
we conclude that all photometric bands show variations. 
Although we can not decide on the basis of Fig. ^whether 
these variations are correlated with the 218 day period in 
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Fig. 2. a) 0-statistics on the 157 data points of the Johnson Vj-band. At least five sub-harmonics can be identified (solid lines). 
On both sides of the main frequency the beat with the year is present (dashed line), b) 0-statistics on the grid of observed 
dates with a synthetic sinusoidal light curve best representing the observations, c) the same as "b" , but on a daily grid (see text 
for explanation). 



the Vj-band photometry it is suggestive that the 218 day 
period is also present in the other photometric bands. 

A powerful tool to study the photometric variations are 
the photometric indices. Both in the Geneva and in the 
Stromgren system there are indices for the effective tem- 
perature Tcff, the Balmer discontinuity, and the blocking 
by absorption lines. 



Bakker et al. (1996) have shown that the photometry 
is consistent with a late B-type supergiant (e.g., B8-9I-II, 
Teff = 12000 ± 1000 K, logg = 3.0 ±1.0, and E{B -V) = 
0.49 ± 0.05). To calculate dereddened photometry we will 
adopt the reddening of E{B - V) ^ 0.49 ± 0.05. Using 
Eq. ^ the Stromgren photometry can be dereddened. No 
dereddening has been applied to the Geneva photometry. 



E{b - y) = 0.70 X E{B -V)= 0.343 ± 0.04 



(6) 



The effective temperature indicators are (6 — y)o 
(Eq. 7) and {B2 - VI) (Eq. 8). 



{b - y)o - (& - y) - E{b - y) 
(B2 - VI) =B2-V1 



(7) 
(8) 



In Fig. ^ both indices are folded with the 218 day period 
and no clear phase dependence is observed. From this we 
conclude that the 218 day period is not related to a peri- 
odicity in Teff. 

The photometric indicators for the Balmer discontinu- 
ity are c? (Eq. 9) and d (Eq. 10). 



c1=^{u~v)-{v-b)- 0.20 X E{b - y) 
d={U -Bl)- 1.430 X {Bl - B2) 



(9) 
(10) 



Fig. H clearly shows the presence of a 218 day period in 
those two indices. The amplitude of the variations in d 
and c\ are about 0.18 mag (a 35(7 detection), i.e. a factor 
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Fig. 3. Top two panels the light curve as a function of phase (left with 11 — 218 and right with 11 — 436 days). Lower panel the 
hght curve as a function of Julian date. In all three cases a sine function is over plotted which does best represent the data for 
the period of 218.0 days (this means two sines per period for 436 days). Geneva data is represented by a dot, and LTPV data 
by a circle. 



eleven larger than the amplitude of 0.016 mag found in the 
Vj-band photometry. The fact that the Geneva d and the 
Stromgren indices are not exactly in phase is likely due 
to the fact that they are not observed simultaneously and 
uses different photometric bands. Furthermore we should 
note again that we did not correct for the long-term pho- 
tometric variability and that this will surely causes an 
additional scatter in the phase light curve. A maximum 
Balmer discontinuity occurs at phase 0phot = 0.0 ± 0.1, 
which is clearly out of phase with the Vj-hand: it is 
d4> — —0.25 ±0.10 behind. The minimum Balmer disconti- 
nuity occurs at phase (/)phot = 0.5 ±0.1, exactly half a pe- 
riod shifted from the maximum Balmer discontinuity. Sup- 
pose that the large standard deviation of the temperature 
indicator (6 — y)Q is due to variations in the effective tem- 
perature of the star. Then the mean relation between Cj* 
and {b — y)o, {b—y)o = —0.116 + 0.097 x c^, predicts a stan- 
dard deviation of the c? index of cr = 0.0008. The observed 
standard deviation is a factor 50 larger and (again) tells us 



that it is not the temperature but the (wind)density that 
causes the variability of the Balmer discontinuity. From 



the work by Bakker et al. (1996) we know that at wave- 
lengths shorter than the Balmer jump (3647 A) the wind 
spectrum dominates and is optically thick in the Balmer 
continuum: the UV energy distribution resembles an A5I 
star. The periodicity in can now be attributed to a 
periodicity in the density of wind material (not the pho- 
tosphere) . 

Photometric variations due to obscuration by wind 
material is also observed in the case of the well studied ob- 
ject Pleione (Luthardt and Menchenkova [1994 ) . The spec- 
tral type of Pleione of B8V (M=3.5 Mq) is comparable 
to the spectral type of HD 101584: B(e)8-9II. In contrast 
to HD 101584 with t;sinz « 50 ± 10 km s^^ (assuming 
that we don't see it pole on), Pleione is a rapid rotator 
(wsinz — 300 km s^^ ) which experiences discrete mass- 
loss episodes. A shell of material is expelled and moves 
away from the star and the spectrum of a normal B8V 
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Fig. 5. Phase light curve of photometric indices. The phase has been calculated by folding the data over a 218 day period and 
taking JD = 2445044.576 as zero point. The three top panel are Geneva indices, while the lower three panels are Stromgren 
indices (left panels: effective temperature; middle panels: Balmer discontinuity; right panels: line blanketing). All sbc abscissae 
are scaled to the same magnitude range of 0.30 mag. 



star changes to a spectrum which is almost completely 
dominated by wind lines which are partially in emission. 
In the case of Pleione the recurrence time scale is of the 
order of several decades, while the recurrence time scale 
for HD 101584 is much shorter (218 days) and the star is 
almost continuously in a shell stage. 

Recent theoretical work on th e shell spectrum after a 
nova outburst (Beck et al. 1995 ) has shown that in the 
case of a dense shell of material around a 10.000 K star, 
the Stromgren sphere of e.g., hydrogen is one order of 
magnitude smaller than in the case of a black-body. This 
leads to a depression of the UV continuum. The observed 
photometric variations of HD 101584 can be understood in 
a similar way as a period depression of the UV continuum 
due to variations in the amount of gas in the line-of-sight. 

The photometric indices for blocking by absorption 
Hues are (Eq. 11) and g (Eq. 12). 



(u - 6) - (6 - y) + 0.18 X Eib - y) 
g^{Bl- B2) - 1.357 x {VI - G) 



(11) 
(12) 



Fig. ^ shows that no clear phase dependence of these 
indices is observed. The presence and strength of spectral 
features does not vary significantly, which indicates that 
the periodicity found can be attributed to a periodicity in 
the strength of the UV continuum level. The number of 
wind lines in the optical spectrum of HD 101584 decreases 
to the red. For wavelengths longer than 5000 A there are 
very few of these lines left and it unlikely that these lines 
produce the periodicity in the photometry. Furthermore 
from Fig. ^ we note that although the Balmer continuum 
shows the strongest variability, there are indications that 
even the Paschen continuum varies in strength. 
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Table 4. Statistics on photometry. 



band 


Ac 


S\ 


M 


Cm 


0'M(A/B) 


U 


3458 


170 


7.779 


0.0710 


0.0074 


u 


3500 


300 


8.606 


0.0684 


0.0348 


Bl 


4022 


171 


7.306 


0.0392 


0.0054 


V 


4110 


190 


7.560 


0.0321 


0.0172 


B 


4248 


283 


6.415 


0.0343 


0.0048 


B2 


4480 


164 


7.890 


0.0309 


0.0052 


f, 


4670 


180 


7.262 


n 0209 


0.0143 


VI 


5408 


202 


7.652 


0.0250 


0.0043 


V 


5470 


230 


6.939 


0.0192 


0.0075 


Vj 


5500 


800 


6.934 


0.0228 




V 


5508 


298 


6.930 


0.0249 


0.0036 


G 


5814 


206 


8.034 


0.0243 


0.0050 


B2 - VI 






0.237 


0.0121 


0.0033 


{b-y)o 






-0.039 


0.0077 


0.0120 


d 






1.308 


0.0345 


0.0072 








0.677 


0.0415 


0.0390 
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-0.066 


0.0091 


0.0050 


ml 






0.037 


0.0157 


0.0150 



4-. 4- Conclusions 

We found a long-term variability (possible period) with a 
typical time scale of 4.7 years. The Balmer (and Paschen) 
discontinuity varies periodically with 218 days and that 
the effective temperature and blocking by absorption lines 
do not vary. This periodicity is therefore attributed to 
changes in the continuum level, and especially in the 
Balmer continuum (UV continuum). As the UV contin- 
uum is much lower than expected for a B8 supergiant we 
argue that the UV continuum is depressed due to the many 
absorption lines from a dense stellar wind and that the 
amount of depression of the UV continuum is period. 

5. Doppler velocities 

5.1. The 218 day period of the Doppler velocities 

An important question in unraveling the photometric pe- 
riodicity is whether or not the periodicity is also present 
in the Doppler velocities of the photospheric absorption 
lines, and if it is, whether the period is 218 or 436 days. 
The small number of data points (only 15) is too small 
to detect a period using the PDM technique and we have 
folded the observations over 218 and 436 days. 

Fig. ^ shows the two possibilities (218 or 436 days). 
On the basis of this figure we are not able to distinguish 
between the two periods, but as we have argued that the 
photometric period is 218 days we will adopt a 218 days 
period in the radial velocities. A sine fit with a period of 
218 days yields an amplitude of 3.0 km s^^ with a resid- 
ual dispersion of 3.4 km s~^ . The fit is not in phase with 
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Fig. 4. Standard deviation of the photometric bands (ctm) as 
a function of central wavelength Ac of the band. Clearly the 
amplitude of the variations increases to shorter wavelengths 
and are in excess of the error. The solid symbols represent the 
observations of HD 101584 (dots: Geneva, squares: Stromgren, 
triangles: Vj-band), and the open symbols the observations of 
the comparison stars crA/(A/B)- 

Table 5. Urad: parameters of best fitted sine function. 
Parameter Best fitted value 



Period U 218.0 ± 0.7 days 

I'System 56.3 ± 1.0 km s~^ 

Amplitude Ki 3.0 ± 3.4 km s~^ 

JD[(t>,^d = 0) 2444990.08 

(/>rad is defined when Uiad = ^system and d«rad/rft > 

c\ but leads it by 0.25 in phase. It is however very doubt- 
ful whether a sine fit is allowed as the average velocity of 
56.3 ± 1.0 km s~^ is 6.0 km s~^ red-shifted with respect 
to the system velocity. Most likely due to the small num- 
ber statistics the observed average radial velocity is off by 
about 6.0 km s~^ . As only one point is below the system 
velocity a sine fit with an offset of the system velocity does 
not not give a reasonable fit to the data. An ad hoc in- 
terpretation to the 6 km s~^ offset could be that we don't 
see direct light from the star, but only reflected light on 
a receding "mirror" of material (e.g., wind material). As 
this is indeed an ad hoc interpretation we will not pursue 
it. 

We have investigated if the periodicity in the Doppler 
velocities could be attributed to changes in the absorption 
line profiles. Fig. shows the line profile of the He I lines 
at six different phases and Fig. || shows the line profiles of 
the two C II lines at 13 different phases. The dashed lines 
is the system velocity of 50.3 ± 2.0 km s~^ with on both 
sides a line at wsin j w 50 ±10 km s^^ to show the variable 
asymmetry of the line profile. The red wing is steeper than 
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Fig. 6. Radial velocity curve of the high-excitation photospheric absorption lines. The two left panels are made by folding with 
the photometric period (218 days), while the right two panels were folded with twice the photometric period (436 days). The 
solid sine is the best fit through the data points (e — 0) and the dashed lines are highly eccentric orbits fitted to the data 
points(left e — 0.7 and right e = 0.2) as a possible explanation for the 6 km s~ velocity offset. The lower two panels give the 
radial velocity curve versus Julian date for 218 and 436 days respectively. C II is represented by a filled circle and He I by an 
open circle. 



the blue wing, and there is clearly a shift in the position of 
the minimum of the profile. This could cause a systematic 
error in the determination of the Doppler velocities. 

5.2. Stellar rotation? 

An estimate of v sin i can be obtained from the full width 
of the absorption line profile of CII (Fig. ||). The full 
width of the He I line is not a good estimate because of 
Stark broadening. From the CII lines we find usini « 
50 ± 10 km s~^ . If 218 days is the stars rotational period 
and the system is almost edge-on (sini — 1) then we find 
that the radius of the star is 214 ± 50 R0. This is one 
order of magnitude too large for a post-AGB star, but is 
about the size expected for a massive 40 M© supergiant. 
Alternative if we adopt a stellar radius (i?, = 20 R©) for 
a post-AGB star, and vsini = 50 ± 10 km s^^ we find a 



rotational period of Prot = 20 days. We can thus conclude 
that 218 days is not the rotational period for a post-AGB 
star (it could be the rotational period if HD 101584 is a 
massive B-supergiant). 

5.3. Pulsation? 

With an amplitude of 6 km s~^ and a 218 day period the 
increase of radius in one pulsation period is 51 R©. This is 
larger than the radius of a post-AGB star (i?* — 20 R©). 
Such an enormous increase in radius would be observed 
and strong variations in the photometric indices which 
are a measure of the effective temperature would be ob- 
servable. As this is not the case we can discard pulsation 
as the origin of the photometric variations. 



Eric J. Bakker et al.: The 218 day period of the pecuhar late B-type star HD 101584 



11 



D 
>> 
D 
X3 

O 
CO 

C\l 



-e- 

+ 



CD 




150 



100 



-50 



AV, 



systen 





[km 



50 



100 



150 



Fig. 7. He I (5876 A) lines profiles as function of phase 
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Fig. 8. C II lines profiles as function of phase ^rad- Each spec- 
trum is offset with 0rad- This figure demonstrated the small 
changes in the position of the center of the profile. The solid 
line is CII (6578.03 A) and the dashed line CII (6582.85 A). 



5.4. Binarity? 



If the 218 day period in the Doppler velocities is attributed 
to orbital motion, the mass function of the system is 

/(to) = 4.9 X 10-^ (1 - e^f^^ Mq. From Bakker et al. 
( 1996 ) we take the system edge-on, sini = 1, and for sim- 
plicity we will adopt e = 0. 



Table 6. Possible 
/(m) = 4.9 X 10~^ M©. 



orbital 



parameters 



for 



Ml = 0.6 Mq 
Ri = 20 Rq 



Ml = 8 Mq 

Ri = 24 Rq 



Ml = 40 Mq 
Ri = 160 Rq 



M2 0.14 Mq 

asini 6.9 R* 
K2 26 km s~ 



0.72 Mq 
13 R* 
67 km s^ 



2.1 Mq 
3.3 R, 

116 km s^i 
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Table | gives the orbital parameters for a 0.6 Mq post- 
AGB star, a 8 M© (Pop. I) and a 40 M© star (Pop I). 
The orbital separation, asinz, is in all cases less than 13 
R*, and all systems are close binaries. In all cases the 
secondary is a low-mass star, with a mass less than the 
primary (the observed star). Stars which are not on the 
main sequence evolve faster than on the main sequence. 
This means that the secondary is most likely either the end 
product of stellar evolution, or is on the main sequence. 
The secondary of the 0.6 Mq primary is a low-mass main 
sequence object as 0.14 Mq white dwarfs do not exist. The 
secondary of the 8 Mq and 40 Mq primary cannot be a 
white dwarf as the progenitor was a high-mass star which 
evolves to a neutron star or a black hole which have a 
mass in excess of 1.4 Mq. The presence of a black hole or 
neutron star is not likely as no high energetic processes 
have been reported in the literature. This secondary is 
thus also on the main sequence. 

A possible solution to the low-mass of the secondary 
and the velocity offset of 6 km s~^ is if the star is in a 
highly eccentric orbit. This will increase the amplitude, 
Kl, the mass-function, and therefore the mass of the sec- 
ondary. Two examples are plotted in the top two panels of 
Fig. ^. In that case we find post-AGB secondaries masses 
in the range between 0.1 and 1.0 Mq. 

It is evident that if the variations in Doppler veloci- 
ties are due to binarity, the mass function poses strong 
constraints on the mass of the secondary and the eccen- 
tricity of the orbit. The most likely evolutionary status, 
the post-AGB phase, yields a mass of the secondary in 
the range of 0.14 Mq for a circular orbit to 1.0 Mq in a 
(highly) eccentric orbit. 

5.5. Conclusion 

Pulsation and rotation are not able to explain both the 
photometric and radial velocity variations. If the observed 
star is in a binary system with an eccentric orbit the radial 
velocity variations could be due to binary motion which 
gives a secondary mass between 0.14 and 1.0 Mq. 

6. Discussion 

We can summarize the results in the following way: 

1. The photometry reveals a long-term variability on a 
time scale of about 1700 days (4.7 years). On the basis 
of the available data we argue that it is likely periodic. 

2. The Johnson Vj-band, derived from the Geneva V and 
Stromgren y-band, shows a highly significant period of 
n = 218.0 ±0.7 days. 

3. The period is most clearly present in the photometric 
indices which are a measure of the Balmer jump. 

4. The effective temperature indicators do not vary with 
phase. 

5. The line blanketing indicators do not vary with phase. 



6. The photometric variations are therefore attributed to 
variations in the (UV and optical) continuum level. 

7. A change in the strength of the Balmer jump can occur 
at a part of the stellar surface with a lower effective 
gravity, leading to a higher mass-loss rate and a flatter 
velocity law of the wind. 

8. The Doppler velocities of high-excitation photospheric 
absorption lines are likely periodic with 218 days, and 
attributed to orbital motion in an highly eccentric or- 
bit. 

9. The Doppler velocities are not in phase with the 
Balmer jump. 

We have discussed stellar rotation, binarity and pulsa- 
tion as possible origin of the 218 day photometric period 
and found that only binarity can explain the observations, 
(a) HD 101584 is a post-AGB star and the variations are 
due to a companion with a period of 218 days. The changes 
in the velocity law and mass-loss rate are likely triggered 
by binary interacting with the secondary. The asymmetric 
photospheric lines profiles might cause a systematic effect 
in the determination of the Doppler velocities, and can 
(partly) explain the 6 km s'^ offset, (b) If HD 101584 is 
a massive object than the stellar rotational period is the 
most acceptable origin. 

In both scenario the same idea works: a part on the 
stellar surface experiences a flatter velocity law and higher 
mass-loss rate than average. This results in a periodic in- 
crease and decrease of the (UV) continuum level. The pe- 
riodicity in Vj-band is also due to the increase in the wind 
density, but as the effect of the amount of neutral hydro- 
gen strongly decrease to longer wavelength the amplitude 
is quite small. 

A alternative explanation for the variations in the 
Balmer discontinuity might be the presence of a second 
continuum source in a binary system (e.g., an accretion 
disk around the secondary). The second continuum source 
is periodically observed or obscured and this leads to the 
observed changes in the photometry. 

From this work we see that the ultraviolet spectral 
region is of crucial importance in understanding this ob- 
ject as the stellar wind is so dominant. We therefore have 
started a long-term lUE program to monitor HD 101584 
every 40 days at high- and low-resolution. Here we predict 
that the UV continuum will rise and fall in phase with the 
Ci and d index. 
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